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Van der Pod's method (Rheol. Acta 1, 198 (1958)) for calculating the shear modulus of a par- 
ticulate composite agrees well with experimental data, but its validity has been questioned, and it 
was applicable only to composites in which the matrix material is incompressible. These limitations 
are removed in this paper in which an error in the original derivation is corrected, and the method 
generalized to apply to any matrix material. Calculations using the corrected theory show that despite 
the error, a table of shear modulus values published with the original theory is sufficiently correct for 
most practical purposes. Applicability of the generalized method to the large class of composites 
having compressible matrices is discussed. Shear moduli calculated by the corrected and extended 
method are compared with corresponding values calculated by other methods currently used. 
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1 . Introduction 

Consider an idealized composite material consisting 
of small spheres imbedded in a matrix. The spheres 
are of approximately the same size, are firmly attached 
to the matrix, and are uniformly distributed so that the 
composite material is macroscopically homogeneous 
and isotropic. An important question then arises — that 
of calculating the two mechanical moduli characterizing 
this material in terms of the moduli characterizing 
the filler spheres and the matrix. 

Several solutions have been given to this problem. 
The best known ones are those of Hashin and Shtrik- 
man [l], 1 Kerner [2], and Budiansky [3]. Hashin and 
Shtrikman's formulae, which provide upper and lower 
bounds for the shear and bulk moduli of the composite 
as a function of filler content, are: for the highest 
lower bound (HLB), 
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for the lowest upper bound (LUB), 
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K m -Kf 3K f +4G f 

In these formulae G represents the shear modulus, 
K bulk modulus, and <p volume fraction, the subscripts 
/ and m refer to the filler and matrix respectively. 

Kerner's formulae, although derived by a different 
method, can be shown by algebraic manipulation to 
be the same as those for Hashin and Shtrikman's 
highest lower bound (1, 2). 

Budiansky provides the pair of coupled implicit 
equations: 
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1 Figures in brackets indicate the literature references at the end of this paper. 
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Young's modulus E and Poisson's ratio v can be 
calculated from the well known expressions: 



Let the cube of homogeneous material (first cube) 
be subjected to the deviatoric stress system 



9KG 
3K+G 

SK-2G 
6K + 2G' 



(9) 



(10) 



<J x —(Ty — —T, (T Z = 2T, Tj:y=Ty Z =T zx =0 
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In van der Pod's solution [4], for which excellent 
agreement with experimental data is claimed [4, 5], 
the formula for K can be shown to be the same as eq 
(2), but the calculation for G is too complicated to be 
expressed conveniently in closed form. Van der Poel 
has treated the special case in which the matrix of 
the composite has a Poisson's ratio ^ = 0.5. The value 
of G is obtained by interpolation from a table of values 
he has provided. It would be useful to have a formula 
applicable when the matrix has any value for Poisson's 
ratio. This was the purpose of the present work. 
During its progress an error was discovered in van 
der Poel's solution. The corrected form is given along 
with details of the derivation. 

2. Derivation 

The derivation is based on a method developed by 
Frohlich and Sack [6]. 

To simplify calculations consider a unit length 
chosen such that on the average there is one sphere 
of filler in each spherical volume 47r/3 of unit radius. 
If a is the radius of the filler sphere, then by defini- 
tion of the unit radius, a 3 is equal to the volume frac- 
tion <P f of the filler. 

Consider a large cube of homogeneous material 
with edges parallel to the axes of a cartesian coordinate 
system having its origin at the center. Consider also 
a second cube in which the region at the center has 
the following structure: A filler sphere of radius a is 
located at the origin. Surrounding this out to a radius 
unity is a sphere of matrix material which in turn is 
imbedded in the homogeneous material. The mechan- 
ical properties of this homogeneous material are 
assumed to be the same as the average macroscopic 
properties that are sought for the composite material 
considered here. 

If the same boundary stresses are applied to each 
of the cubes, it is assumed that the displacements in 
the two cubes are the same at a distance r = R > 1 
except for terms of a high order in 1/r. In the second 
cube the influence of the structure within the radius 
r= 1 on the displacement at r=R is expected to be of 
the same order as the ratio of the two volumes, i.e., 
~ l//? 3 . Thus it will be required that 

Displacement at r=R in second cube 
Displacement at r=R in first cube 



l +l 



const 
R" 



(11) 



where the as represent tensile and the rs shear 
stresses. Assuming Hooke's law, it is easily seen that 
this stress system gives rise to the following displace- 
ments u x , u y , u z within the material; 



TQ L +v)x = Tx 

E " 2G 
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In spherical polar coordinates the components of 
stress and displacement are obtained by transforma- 
tion as 



O- r =2TP 2 (Cos0), Tr0=TPL(Cos0), Trp = 



(14) 



Tr Tr 

u r = -77 P 2 (Cos0), ue=-r^, 



P! z (CosO), i^ = (15) 



where P 2 (Cos0) is the second Legendre polynomial 
and P^Cosfl) is its derivative with respect to 0. 



P 2 (CosO)=i(3Cos 2 6-l), 
P.; (Cos 6) =-3 Cos 6 Sin 6. 



(16) 



In order to find the displacements in the cube con- 
taining the special spherical structure at the center 
(second cube), it is necessary to find appropriate 
solutions to Lame's equations expressed in spherical 
polar coordinates. These equations are tabulated by 
Love [7]. Because of symmetry of the stress system 
about the z axis, the displacements u r and ue are 
independent of the azimuthal angle <p, and the dis- 
placement u<p = 0. Lame's equations thus simplify to 



(A + 2G) Sin0^ + 2G-^- (rco*Sin0)=O 
36 dr 

(A + 2G)r Sin 6 ^~ ~ 2G 4h («*> Sin 6) = 

or ou 



(17) 



where according to Love [8] the dilatation A is given by 
1 



= ^ti<^ Sin *> + l><™« s H 



and the rotation a><p by 



0)<p 
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2r\dr 
Lame's constant X is given by 
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The most general solution having the same P 2 
symmetry as the displacement system (15) is obtained 
by setting u r = r$P 2 (Cos 0), uo = arPP 2 ' (Cos 6) in eqs 
(17) to (19), and substituting for k using eq (20) to obtain 
the set of linear homogeneous equations in a: 

[(l-2i>) (j3+l)j3-12(l-i>)]a+ [/3+4(l-*/)] = 

(21) 
[-3£ + 9-12*/]a+ [(1-i>)03+1)/3-5 + 8i>] = O. 

In order that these equations be compatible, it is 
necessary that 

+ 4(1 -v) 



{\-2v){£+ 1)^-12(1 -v) 

(1 -v)(p+ l))8-5-h8i/ 



(22) 



— 3/3 -h 9 — 12^ 

or, after simplification 

0S-l)(/3 + 2)08-3)03 + 4)=O (22) 

The value of a corresponding to each root of eq (22) 
is found by substituting into eq (21). The solutions for 
u r and ue thus become. 2 

u r = [Ar+Br~* + Cr* + Dr~ 2 ]P 2 (Cos 6) (23) 

J P' 2 (Cos 0) (24) 



uo~- 
where 



\Ar-\Br* + bCr* + cDr 2 



7-4* ; 
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» c r 



\-2v 
5-4*/ 
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and A, B, C, D are constants. 

The strains e r r and e r e are found from the relations 



fhere 
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Within the filler sphere r ^ a the requirement that 
the displacement be zero at the origin necessitates that 
B/= D/= 0, where the subscript / refers to the filler 
medium. The solution applicable in the region a ^ r ^ 1 
involves the constants A m , B m , C w , D m where the sub- 
script m refers to the matrix medium. Within the 
region 1 ^ r ^ R requirement (11) necessitates that 
A = T/G, C = D = 0. A must be regarded as an unknown 
constant as it involves the shear modulus G, the quan- 
tity which is to be calculated. 

At the boundaries r=a and r= 1 the radial and 
tangential components of the displacement and the 
stress must be continuous. This gives rise to the 
following set of linear homogeneous equations in the 
unknowns Af, Cf, A m , B m , C m , D m , A, B. 



MA f ~ Ma 2 C f -A,„ + 4a-*B m + |a 2 C„, - k 3m a^D m = 



MA f + Mk 2f a 2 C f -A„,-^a^B„ 



k>mCl C m fi4mtt ' Dm — U 



<iA/-\- d*Cf— aA m — a~ 4 B m — a 3 C m — a *D„, = 



- aA f + b f a 3 C f — - aA m +- a- 4 B n 



- b m a 3 C m — c m a~ 2 D m = 0. 
(31) 



d d Uo , 1 d 

--u r ,e r e = -ue-- + - Ve u r (26 ) 



A m - 4B m -- C,„ + W>* + 4W£ - NA = 



and the stresses necessary for the solution are then 
calculated using 



26V 



+ 2Ge r r? Tr6—Ge 
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(27) 



1 \-2v 
to obtain 2 
o- r = 2G[A-4Br-* + k ] Cr 2 + k,Dr-z]P 2 (Cosd) (28) 



T r — G 



^+|^r- 5 + A 2 Cr 2 + 



Dr" 3 1 P; 



(Cos 0) (29) 



2 For an incompressible material (u — 0.5) van der Poel's solution agrees with the solution 
obtained here, but when v / 0.5 the exponents in the terms having coefficients C and D in 
van der Poel's solution are erroneously given as functions of v. This error has previously 
been noted by Hashin [9]. 



A m + - B m + k 2m Cm + k Am D m --NB-NA = Q 
A m + B m + C m + D m -B-A = 



-A„, — - B m + b m Cm + c,„D m + -B —-A = 



yhere 



M = G f IG m ,N=GIG m . 



(32) 



In order that the solution be nontrivial the deter- 
minant of the coefficients must be set equal to zero. 
Solving the resulting equation will evaluate the shear 
modulus of the material G, which is the quantity 
desired. 
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After simplifications the determinantal equation 
becomes 3 



2M 
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(34) When the volume fraction <p ■ 


= a 3 of the 



(33) 



(35) 



where X 

Equation (33) can be expanded into the form 

AX 2 + BX + C = 

where 

A = 2(M-l)[M(7 + 5v f ) 

+ 4(7-101//)] [4(7-101/,,,) (4-5z, m )(l-a 3 )(l-a 7 ) 

-63a 3 (l-a 2 ) 2 ] 

+ 280(M-l)(7-10^)(l-i/ m )(4-5^/ m )(l-a 3 )a 7 

+ 60[M(7 + 5*//) 

+ 4(7-10^)](l-^)(7-10^)(l-a 7 ) 

+ 2100(7- 10i/,) (l-^) 2 a 7 (36) 

S=70(M-l)[M(7 + 5*//) 

+ 4(7-10^)](l-i/ m )[(4-5i/,„)(l-a 3 ) 

-6a 3 (l-a 2 )] 

+ 525[M(7 + 5i//) + 4(7-10i//)](l-j/ m )2 

+ 60(M-l)[M(7 + 5*//) 

+ 4(7-10i/ / )](l-i/,„)[10i/ m a 3 (l-a 7 ) 

-7a 5 (l-a 5 )] 

-2100(M- 1) (7- lOi//) (1 - */,„) 2 a 10 



C=-525(M-l)[M(7 + 5i//) 

+ 4(7-10*/ / )](l-j/ m ) 2 a 3 . 



(37) 



(38) 



One of the roots of eq (35) is negative and is extra- 
neous. The positive root when substituted into eq 
(34) provides the value of the shear modulus G. 



3 The determinant given by van der Poel agrees with the determinant, eq (33), when 
v m = 0.5 except for elements in the second column. These elements in van der Poel's 
solution are more complicated functions of v j . 



is small, terms of higher order than a 3 can be neglected. 
The expressions for A, B, and C become much shorter, 
and a further simplification results because the term 
[M(7 + 5j//)+4(7 — lOvf)] can be factored from each 
of these quantities. The solution of eq (35) for this case 
yields the following approximate formula for G: 



-=1 + 



15(1— V m )(Gf— G m )(f 

(8-\0p m )G f +(7-5v m )G n 
3. Discussion 



(39) 



It can be shown by algebraic manipulation that the 
formulas (39), (1), (5) and (3) for G can be displayed in 
the form 

G=G m + (G f -G m )$= (l-®)G m + (!>G f (40) 

which, when <& = (p, is the well-known rule of mixtures. 
For the approximate van der Poel equation (39) 



, _ (8 — 10^ m )G m + (7 — 5v m )G m 
(8-l0p m )G f + (7 -Sv^Gm*' 

For Kerner, or Hashin and Shtrikman HLB, (1) 

(8 - \0v m )G + (7 - 5u m )G m 

(8-10^)G / +(7-5i/ m )G m ^- 



(41) 



(42) 



For Budiansky (5) 
(8- 



<D = - 



10v)G+(7-5v)G 



<P- 



(8-10v)G f +(7-5v)G 

For Hashin and Shtrikman LUB (3) 

(8-lQ Vf )G+(7-5 Vf )G f 
(8-Khy)C,+ (7- SvJGf*- 



(43) 



(44) 
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This display shows that these formulas are all of the 
same type, and can be interchanged into each other 
in a certain predictable way. 

When the volume fraction of filler <p is very small so 
that G ~ G m , and v ~ v m , the values of <t>/^> computed 
by eqs (41) to (43) are essentially equal. Thus plots of 
<I> = (G — G m ) I (Gf— Gm) versus tp for these relations all 
have the same initial slope. Plots of GjG m for these 
relations also have the same initial slope. Thus the 
initial increases of GjG m with increasing <p, as calcu- 
lated by the methods of Hashin and Shtrikman HLB, 
Kerner, van der Poel, and Budiansky, are the same. 

The slope G/G m computed from the approximate 
van der Poel formula (39) remains constant as <p 
increases, but the slope computed by the accurate 
method (35) is nonlinear because of the complicated 
dependence on various powers of a the filler sphere 
radius. The slopes computed using the formula of 
Kerner, or Hashin and Shtrikman HLB, (1) and 
Budiansky (5, 6) are also nonlinear, and increase with 
increasing^. 

When v m = 0.5 and GjjG m > 1 as in the case of a gum 
rubber matrix loaded with rigid filler particles, the 
shear modulus for small values of <p should obey the 
equation 



Table 1. G/G m for various values of v { and G f /G r 
according to Budiansky's formulae (5, 6) 



calculated 



G=G m (l + 2.S<p). 



(45) 



This equation is analogous to an equation developed 
by Einstein [10] for the viscosity of a suspension of 
rigid spheres in a fluid medium, and has been used to 
predict reinforcement effects in rubber [11]. Under 
similar conditions the van der Poel formulae developed 
here (35), (39) and the formulae of Kerner, or Hashin 
and Shtrikman HLB, (1) and Budiansky (5, 6) reduce 
to eq (45). In fact the approximate formula (39) may be 
considered as a generalized form of eq (45) in which 
Gf and v m are specified in addition to G m . 

In order to calculate G using the formula of Kerner, 
or Hashin and Shtrikman HLB, (1), it is necessary to 
know only one of the moduli characterizing the filler 
material, i.e., Gf. This is in contrast to Budiansky's 
formulae (5, 6), where elimination of K leads to a 
relationship for G involving filler moduli Gf and 
Kf. Similarly for the corrected van der Poel computa- 
tion (35), two filler moduli Gf and Vf are required. 
The corrected van der Poel computation (35) and the 
formulae of Budiansky (5, 6) thus provide some indi- 
cation of the relative effectiveness of the distortional 
(shear modulus) and dilatational (bulk modulus) 
properties of the filler in determining the shear modu- 
lus of the composite. To determine this relative 
effectiveness, computations were made by the two 
methods and the results are presented in tables 1 
and 2. 

In table 1 values oiGjG m calculated from Budiansky's 
formulae (5, 6) are given, corresponding to various 
values of filler Poisson's ratio Vf and the ratio G//£ m , 
for a matrix Poisson's ratio v m = 0.5. In the calculations 
the ratio of the shear or rigidity modulus of the filler 
to that of the matrix GjjG m ranges in value from 10 to 



Vf 



GjGm 

corre- 
sponding 

to 
G f /G m = 10 



100 



1,000 



10,000 



100,000 







<p = 0.25, 


i/ w = 0.50 






0.10 
.30 
.50 


1.7759 
1.7861 
1.7958 


2.4647 
2.4726 
2.4787 


2.6428 
2.6440 
2.6449 


2.6642 
2.6644 
2.6644 


2.6664 
2.6664 
2.6664 






V? = 0.50, 


i/ m = 0.50 






0.10 
.30 
.50 


3.4855 
3.5647 
3.6565 


17.871 
19.003 
20.143 


146.56 
158.64 
170.74 


1426.2 
1548.6 
1670.8 


14221. 
15448. 
16671. 



Table 2. G/G m for various values of v { and G f /G m , calculated 
according to the corrected van der Poel equation (35) 



Vf 


G/Gm 

corre- 
sponding 

to 
G//G m =10 


100 


1,000 


10,000 


100,000 






<p = 0.25, 


I'm = 0.50 






0.10 
.30 
.50 


1.6842 
1 .6843 
1.6845 


2.0025 
2.0026 
2.0027 


2.0490 
2.0491 
2.0491 


2.0540 
2.0540 
2.0540 


2.0545 
2.0545 
2.0545 






<p = 0.50, 


i/ w = 0.50 






0.10 
.30 
.50 


3.1083 
3.1235 
3.1393 


6.4334 
6.4888 
6.5356 


7.8687 
7.8810 
7.8909 


8.0738 
8.0751 
8.0762 


8.0952 
8.0953 
8.0954 



10 5 , and from eq (10), the range of Vf from 0.1 to 0.5 
corresponds to a range of the bulk modulus Kf from 
(11/12)G/ to oo. The tabulated values show that the 
relative shear modulus GjGm (Budiansky) can be 
strongly dependent upon the rigidity of the filler relative 
to the matrix G f /G m , but is only slightly sensitive to the 
value of the filler Poisson's ratio Vf or the bulk modulus 
K f . 

Table 2 presents similar results calculated using the 
corrected van der Poel formula (35). Here again the 
values of the relative modulus GjG m are seen to be 
insensitive to the value of the filler Poisson's ratio. 
Thus although the Budiansky and corrected van der 
Poel formulas make use of the value of filler Poisson's 
ratio in the calculations of G\G m , the values obtained 
are insensitive to it, and in this respect these formulas 
offer little advantage over the formula of Kerner, or 
Hashin and Shtrikman HLB, in which the value of 
filler Poisson's ratio is not needed. 

Table 3 presents values of G/G m for various filler 
volume fractions <p and values of GflG m , computed 
according to the corrected van der Poel equation (35). 
In these calculations Poisson's ratio of the filler spheres 
was ^=0.25, and of the matrix was v m = 0.5. If these 
values are divided by GfjG m the resulting values are 
almost identical with the values of GjGf given in a 
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Table 3. G/G m as a function of<pfor various values of G f /G m , calculated according the corrected van der Poel equation (35) 









I'm = 0.50 


, v f = 0.25 










G/Gm 














<p 


corresponding to 
GflG m = z 2 


10 


20 


100 


1,000 


10,000 


100,000 


.00 


1.0000 


1.0000 


1.0000 


1.0000 


1.0000 


1.0000 


1.0000 


.05 


1.0363 


1.1025 


1.1166 


1.1296 


1.1328 


1.1331 


1.1332 


.10 


1.0738 


1.2177 


1.2504 


1.2814 


1.2891 


1.2899 


1.2900 


.15 


1.1126 


1.3499 


1.4087 


1.4665 


1.4813 


1.4828 


1.4829 


.20 


1.1529 


1.5038 


1.6000 


1.6998 


1.7262 


1.7289 


1.7292 


.25 


1.1946 


1.6843 


1.8350 


2.0026 


2.0491 


2.0540 


2.0545 


.30 


1.2379 


1.8962 


2.1264 


2.4066 


2.4901 


2.4991 


2.5000 


.35 


1.2825 


2.1438 


2.4888 


2.9606 


3.1162 


3.1333 


3.1350 


.40 


1.3286 


2.4301 


2.9376 


3.7398 


4.0447 


4.0798 


4.0834 


.45 


1.3761 


2.7559 


3.4852 


4.8580 


5.4924 


5.5713 


5.5794 


.50 


1.4249 


3.1196 


4.1376 


6.4759 


7.8782 


8.0748 


8.0953 


.55 


1.4750 


3.5174 


4.8898 


8.7861 


12.031 


12.578 


12.636 


.60 


1.5264 


3.9448 


5.7280 


11.948 


19.569 


21.236 


21.423 


.65 


1.5791 


4.3986 


6.6336 


15.978 


33.437 


38.824 


39.482 


.70 


1.6332 


4.8794 


7.5946 


20.684 


58.323 


76.204 


78.710 


.75 


1.6889 


5.3937 


8.6140 


25.740 


99.923 


159.63 


170.20 


.80 


1.7463 


5.9570 


9.7226 


30.888 


160.66 


356.81 


408.82 



similar table by van der Poel [4]. The elements in the 
second column of the determinantal equation evaluated 
by van der Poel involved different functions of vj 
than those given in eq (33) of this paper. Evidently van 
der Poel's results agree with those given in table 3 
because of the insensitivity of the method of calculation 
to values of vj. 

The way in which results of the corrected van der 
Poel calculation compare with results calculated by 
the other methods considered here is illustrated in 
figures 1 and 2. In the example chosen for figure 1, 
Poisson's ratio of the matrix is ^ m = 0.4 and the ratio 
of filler to matrix rigidity is G//G m = 30. For those cases 
in which it is used, the value of Poisson's ratio of the 
filler is z^=0.25. The plots of relative shear modulus 
GIG m versus <p represent the expected shear properties 
of a series of particulate composites consisting of 
small glass spheres imbedded in a rigid epoxy matrix. 

In the example depicted in figure 2, v m = 0.5 and 
G f /G m = 70,000. When used, iy=0.25. The plots in 
this case represent the expected shear properties of 
a series of particulate composites consisting of small 
glass spheres imbedded in a lightly vulcanized matrix 
of natural rubber. 

For the epoxy matrix composite, figure 1, the pre- 
dictions of the van der Poel and Kerner, or Hashin 
and Shtrikman HLB, equations are about the same, 
but for the rubber matrix composite, figure 2, the 
predictions of these two formulas are much different. 
The Budiansky calculation predicts the highest 
values of GjG m for both examples. 

From assumptions made in the derivation, one 
should only expect reliable predictions from the van 
der Poel calculation when the volume fraction of the 
filler <p is small. However, agreement of the theory 
with experimental results appears to be good even for 
large values of the volume fraction of filler. According 
to van der Poel [4], calculations for the case 



G f IG m = 100,000, 1^=0.25, i/ m = 0.5 



are in good agreement for values of <p up to 0.60 with 
Eilers' [12] experimental values for the relative viscos- 
ity of suspensions of bitumin particles in water. 

If the values in table 3 corresponding to these van 
der Poel calculations are compared with Eilers' data, 
the same degree of agreement is found. In table 3 
there is only slight difference between the values for 
Gf/G m = 10 4 and GfjG m = 10 5 . Therefore the curve in 
figure 2 of GjGm versus <p, calculated using Gf/G m 
= 70,000 in the corrected van der Poel formula (35), 
lies very close to a similar curve that could be plotted 
from Eilers' data. The curves plotted using Budiansky's 
formulae and the Kerner, or Hashin and Shtrikman 
HLB, equation are quite different. It seems then that 
either eq (35) or van der Poel's original calculation 
provides the best agreement with Eilers' data. 

Schwarzl [5] has studied the mechanical character- 
istics of composites consisting of NaCl crystals in a 
rubbery polyurethane matrix. He reports agreement for 
volume fractions of filler up to <p = 0.50 between his 
results and van der Poel's theory using the parameter 
values Gf/G m = 10 4 , i^/=0.25, and i^ m = 0.50. The curve 
in figure 2 calculated from the corrected van der Poel 
equation provides a good fit to Schwarzl's data, and 
the other two curves do not fit. 

When measurements were made at low temperature 
so that the polyurethane matrix was in the glassy 
state, Schwarzl obtained data that could be fit by 
van der Poel's theory using the parameter values 
Gf/Gm — 8.4, i^/=0.25, and i^ m = 0.5. The value of 
^ m = 0.5 was used because van der Poel's table had 
been calculated for this value only. The value of 
v m that should have been used was probably close to 
^ m = 0.4. This suggests that values of G/G m calculated 
by the van der Poel method are not sensitive to the 
value of the matrix Poisson's ratio i> m , so calculations 
were made to check this. 

In table 4 values of G/G m calculated from the cor- 
rected van der Poel equation (35) are given, corre- 
sponding to various values of matrix Poisson's ratio v m 
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FIGURE 1. Plot of relative shear modulus G/G m versus volume fraction 
of filler for the case G f /G m = .?0, v { = 0.25, v m = 0A0. 

Curve 1, Approximate van der Poel equation (39); 2, equation of Kerner, or Hashin and 
Shtrikman HLB, (1); 3, corrected van der Poel equation (35); 4, Budiansky's formulae 



Table 4. G/G m for various values of v m and G f /G m , calculated 
according to the corrected van der Poel equation (35) 
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FIGURE 2. Plot of relative shear modulus G/G m versus volume fractio 
of filler for the case G f / G m = 70,000, v { = 0.25, v m = 0.50. 

Curve 1, Approximate van der Poel equation (39); 2, equation of Kerner, or Hashin and 
Shtrikman HUB, (1); 3, corrected van der Poel eqi ' 



(5, 6). 



equation (35); 4, Budiansky's formulae 



G/Gm 

corre- 
sponding to 

Gf/Gm = 3 



10 



30 



100 



1 ,000 



10,000 



= 0.25, 



vs 



= 0.25 



0.30 
.40 
.50 



1.2955 
1.3058 
1.3224 



1.5556 
1.5992 
1.6843 



1.6747 
1.7445 
1.8982 



1.7245 
1.8077 
2.0026 



1.7451 
1.8343 
2.0491 



1.7472 
1.8371 
2.0540 



<p = 0.50, ^=0.25 



0.30 
.40 
.50 



1.6928 
1.7156 
1.7542 



2.5562 
2.7189 
3.1197 



3.0870 
3.4400 
4.7709 



3.3470 
3.8359 
6.4759 



3.4626 
4.0231 
7.8782 



3.4747 
4.0432 
8.0748 



and the ratio G f JG m , for a filler Poisson's ratio ^=0.25. 
These calculations show that when G f /G m ~ 10, 
changing v m from 0.5 to 0.4 changes the value of 
G/Gm obtained only by a small amount, which approxi- 
mates the experimental error expected in a measure- 
ment. When GfjG m ^ 30, the change in the calculated 
value of G/G m is significant, so that values of G\G m 
taken from a table calculated assuming ^ m = 0.5 would 
not provide a good fit to experimental data. Thus 



although Schwarzl was able to fit his low temperature 
data using van der Pod's table, the corrected and 
extended van der Poel equation (35) would be needed 
to fit data for a system in which Gj\G m ~ 30, v m = §A 
such as glass spheres imbedded in a rigid epoxy 
matrix. 
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